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Abstract. The Kolmogorov-Sinai (KS) entropy in turbulent diffusion of magnetic 
field lines is analyzed on the basis of a numerical simulation model and theoretical 
investigations. In the parameter range of strongly anisotropic magnetic turbulence the 
KS entropy is shown to deviate considerably from the earlier predicted scaling relations 
[Rev. Mod. Phys. 64, 961 (1992)]. In particular, a slowing down logarithmic behav- 
ior versus the so-called Kubo number R ^ 1 {R = {SB/Bo){^\\/£,±), where SB/Bq 
is the ratio of the rms magnetic fluctuation field to the magnetic field strength, and 
and are the correlation lengths in respective dimensions) is found instead of a 
power-law dependence. These discrepancies are explained from general principles of 
Hamiltonian dynamics. We discuss the implication of Hamiltonian properties in gov- 
erning the paradigmatic "percolation" transport, characterized by i? ^ cx3, associating 
it with the concept of pseudochaos (random non-chaotic dynamics with zero Lyapunov 
exponents). Applications of this study pertain to both fusion and astrophysical plasma 
and by mathematical analogy to problems outside the plasma physics. 
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1. Introduction 

The presence of low-frequency, long-wavelength fluctuations in hot magnetized plasmas 
is often found to alter the plasma confinement properties leading to an anomalously high 
heat and energy transfer across magnetic field lines as compared to purely collisional 
values, a phenomenon known as anomalous or turbulent transport. While the basic 
deteriorating effect of the fluctuations on the plasma confinement is relatively well 
understood, a detailed, microscopic picture of anomalous transport turns out to be 
far than simple. Situations exist in which collisional and fluctuation induced properties 
operate essentially on the same footing, as for instance some regimes of electron outward 
fluxes in solar coronal loops as discussed by Galloway et al. [1]. Well before that, in 
the fusion field, essentially the same mechanism has been proposed for electron heat 
transport due to parallel thermal conductivity in a tokamak with destroyed magnetic 
surfaces |2]. Early investigations of this problem are due to the pioneering work of 
Rosenbluth et al. [3j. In all these models the fluctuations are such as to introduce 
some form of braiding of magnetic field lines [1], with "time" associated with the 
axial (toroidal) dependence. In the braided magnetic fields, the collisional transverse 
diffusion is naturally enhanced as a consequence of the parallel motion providing a 
continuous mapping of the particle's perpendicular random walk. In realistic tokamaks, 
in view of the many different processes actually producing the transport (e.g., Refs. 
El [TJ [HI [9l [To], just to mention some), effects caused by magnetic braiding can, 
however, represent only a fraction of the overall transport of energy and heat and be 
significant within only a limited range of parameters. Nevertheless, addressing the 
magnetic braiding phenomenon in the appropriate setting |2l [3] defines a problem of 
fundamental and practical importance. Not only does it relate to a fascinating set of 
nonlinear dynamics issues [Uj, but also a reliable extrapolation of braiding induced 
transport losses to reactor-relevant plasmas is not at hand. 

In order to accurately assess the collisional transport in the braided magnetic field, 
a few parameters are needed. These parameters are: (i) parallel (^y) and perpendicular 
(^x) fluctuation correlation lengths; (ii) collisional parallel and perpendicular particle 
diffusivities: these quantities are actually known as functions of plasma density and 
temperature in both axial and toroidal geometries [12] ; (iii) fluctuation induced magnetic 
field line diffusion coefficient (-Dm); and (iv) the rate of exponential separation of initially 
close magnetic field lines, h, a parameter often referred to as the Kolmogorov-Sinai (KS) 
entropy [m [IS [H [15] . 

Most previous theoretical and numerical studies have focused on addressing the 
scaling laws for the diffusion coefficient Dm as a function of the fluctuation strength or, 
more conveniently, the Kubo number R = (5i?/i?o)(^||/^±), where 6B/Bq is the ratio 
of the rms magnetic fluctuation field to the background magnetic field. In general, the 
results of those investigations can be summarized by a power law dependence D^, oc W 
(Refs. [ini [171 [HI [IS]), where the scaling is quasilinear-like (7 = 2) for _R <^ 1 and 
percolation-like (7 = 7/10) for R ^ 1 (as distinct from traditionally believed, Bohm- 
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like scaling with 7 = 1 fIU[ f2T\). More recently, the percolation scaling with 7 = 7/10 
was questioned by Milovanov [221 [23], who argued that 7 = 2/3 instead. Support to this 
value can be found in the numerical simulation results of Ref. [211 (25]. Milovanov has 
also addressed [23] the theoretical foundations of the percolation scaling by relating it to 
a description in terms of fractional derivative equations and Hamiltonian pseudochaos 
(random non-chaotic dynamics with zero Lyapunov exponents) [26]. This connection 
to pseudochaotic properties of percolation will play a central role in the theoretical 
investigations below. 

Actually less progress has been made with respect to addressing the h value. Here, 
too, scaling relations have been proposed [161 [27]: h oc R"^ for R <t: I and h oc R^'^ InR 
for R ^ 1. Even so, their direct numerical investigation has been performed only 
recently [28] . A net result of this investigation follows: While the quadratic behavior of 
h in the parameter range of small Kubo numbers R ^ 1 could readily be reproduced, a 
strong deviation from the predicted square root-like scaling h oc R^^^ In R was found for 
R^ 1. Indeed, the observed behavior was markedly slower (Fig. 4 in Ref. [28]), possibly 
slower than a power law, with a tendency to saturation for the largest achievable R (of 
order 10^). This study is aimed to elaborate those investigations and to comprehend 
the observed discrepancies. 

The rest of the paper is organized as follows. A numerical study of the problem 
is presented in Sec. 2, where we detail the numerical techniques and also demonstrate 
the phenomenon of saturation of the KS entropy for the largest achievable R. These 
results are theoretically analyzed in Sec. 3 on the basis of a Hamiltonian approach. 
In particular we discuss the implication of fractal geometric properties in the entropy 
growth slowing down. We summarize our conclusions in Sec. 4. 

2. Numerical study 

We begin with the standard equations for magnetic field lines, with r a three-dimensional 
(3D) position vector, written as 

dT/ds = 6B{r)/\B{r)\, (1) 

where s is the natural parameter, such that / ds is the length along a field line, and B(r) 
is the magnetic field. The latter is set as a sum of a constant homogeneous background 
field Bo = Bqz in a straight-cylinder geometry (here z is a unit vector in the axial 
direction) and a static magnetic perturbation, defined by its Fourier series expansion as 

5B(r) = E E ^^keL'^^ exp z[k ■ r + (2) 

a=l,2 k 

where e^T'' are two polarization unit vectors (these corresponding to the shear Alfven 
and compressional Alfven waves), (j)^'^ are randomly chosen phases, and (5i?k are the 
spectral amplitudes. Basically, this is the same realization of magnetic fluctuations as 
in Refs. [I9l[28]. It is assumed that (5S_k = SB]^ and = so that the field (5B(r) 
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be real. Finally, the fluctuation spectrum is taken in the power law form to be, with 
O = 6(^±,^||) a normalization constant. 

Such power law spectra of varying anisotropy have been found in many plasmas, both 
in space and in the laboratory [291 EOl EU [321 [331 [Ml [35]. In the simulation a = 5/3 
reminding of the Kolmogorov spectrum of fully developed turbulence. It is worth 
to remark that, while the exact entropy value in the simulation may be somewhat 
dependent on a, we expect that the basic physics conclusions regarding the functional 
behavior of h will remain essentially the same. This expectation relies on the results of 
theoretical analysis below, which will lead to (i) generic functional form for the entropy 
in Hamiltonian dynamics, and (ii) comprehension of the features of nonlinear saturation 
of the h function in association with the natural dependency of resonant and percolation 
properties on the Kubo number. 

Note that the numerical model setup in Eqs. ([I]) — ([3]) is essentially 3D. The values 
of the correlation lengths ^y, ^± determine the shape of the constant intensity ellipsoids 
in k-space. Special cases of this correspond to an isotropic turbulence = ^\\); a 
quasi-2D turbulence (^y ^ ^_|_); and the slab turbulence (^y <^ For the present 
analysis, ^\\ > but other cases have been considered in Refs. [T9l [28]. The wave 
numbers are taken from A^^in — + + "-^ ^ ^max; with the wave vectors defined as 
k = (27r/A^min) (nx/^±,ny/^±,nz/^\\^ . To match with the natural physics and computer 
limitations, the power law dependence in Eq. ([3|) is truncated on both sides implying 
finite cutoffs at fc_Lmin = 27r/,^^ and fc^max = (27r/,^^)(A^inax/A^min) in the perpendicular 
plane, and at /c||min = 27r/,^|| and fc||max = (27r/^||)(A^max/A^min) in the axial direction. 
It is convenient to use ^ir/k^^^^ as the unit length in the perpendicular plane. Hence 
the perpendicular correlation length can be expressed as = A^max/A^min- To match 
with the expected properties of heat transport in coronal loops and pertinent model 
parameters of Ref. |28], here we set A^min — 2.23 and A^max = 16, which generates 
some 8200 independent wave modes in Eq. ([2]). We take notice of the fact that the 
KS entropy may depend on the spectral extension in Eq. Q since smaller wavelengths 
favour faster separation of the field lines [371 [3HI [39]; however, for the runs presented 
here, we keep the spectral extension fixed (i.e., N^^^^/N^i^ is not varied). The sizes 
(in respective directions) of the simulation box are set to L^^y = A^min^± = A^max and 
Lz = A^min'Cll = ('C||/^±)^max5 SO that at least two correlation lengths have to be traveled 
before the same magnetic configuration is reproduced. This effectively eliminates the 
periodicity effects [36) . 

In the basic theory [TTl [ID], the KS entropy is obtained as the sum of the positive 
Lyapunov exponents for the unstable motions. The Lyapunov exponents, in their turn, 
can be calculated numerically by implementing the algorithm of Wolf et al [H], in which 
they are obtained from the Gram-Schmidt orthonormalization. In the numerical model, 
6B/Bq is varied from 0.01 to 1, and from 1 to 100; this choice of parameters has 
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Figure 1. The three Lyapunov exponents computed as a function of time for R = 40. 
After an initial transient, the exponents settle to the asymptotic values. 



earlier been discussed for astrophysical applications in Refs. [211 [39l |12]. With these 
settings, the Kubo number R can be varied from 0.01 to 100. These values extrapolate 
from the quasilinear to the percolation regime as discussed in Refs. [2S1[IS]- Also, it was 
shown in Ref . [28] that the value of h depends only on R, and not separately on 6B / Bq 
or ^\\/^±- this work, focusing mainly on the percolation regime, we have restricted 
our investigations to the parameter range R> 1. Overall, from 40 to 100 different initial 
conditions for the field line Eq. ([1]) have been used to compute the Lyapunov exponents 
numerically. The numerical accuracy has been checked by asking that the sum of the 
three Lyapunov exponents be zero, as implied by V ■ B = 0. After an initial transient, 
this was soon achieved with an accuracy of the order of 10~^. We also checked that 
the numerical algorithm has actually reached the asymptotic regime: Figure 1 shows 
the three Lyapunov exponents for i? = 40 as a function of the integration time. It can 
be seen that only one Lyapunov exponent is substantially larger than zero, and that 
the sum of the three exponents equals zero. Also, some small fluctuations are seen in 
the long time regime: this allows to estimate the statistical error of the Kolmogorov 
entropy, which turns out to be less than 0.1 ^J^. The typical results of this investigation 
obtained for ^||/^± = 100 and 1 < i? < 100 are shown in Figure 2. 

There are two main conclusions to be drawn from the analysis. First, for not 
too large the R values, 1 <~ R <~ 30, the h function is well approximated by a 
logarithmic dependence h ~ jjlni? + g, with p and q the numerical fitting parameters. 
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This is shown in detail in the inset of Figure 2, where the fit is made only for the points 
up to i? = 40. Second, with the increasing R the entropy logarithmic growth slows down 
showing signatures of saturation, consistently with our first results in Ref. [28]. This is 
clearly seen in Figure 2, as fitting all the points with the logarithmic dependence shows 
a discrepancy which goes beyond the numerical error bars. When larger R are achieved 
in the simulation, the slowing down effect is more pronounced. These behaviors are in 
marked contrast with the prediction of Ref. [16] that h oc R^^'^ hiR for R^ 1. 

3. Theoretical model 

3.1. Hamiltonian formulation 

In the parameter range of practical importance, the rms magnetic fluctuation is supposed 
not to exceed the background magnetic field Bq. To simplify the analysis, we shall 
assume that 5B/Bq <C 1. Being somewhat restrictive, this condition, however, still 
retains the essential physics of the model. Since we are interested in the regime when 
R^ 1, we assume that ^ 1 with a large margin. Thus, it is required that the 

turbulence be strongly anisotropic and its wave vector spectrum be actually squeezed 
to the perpendicular plane. With this gross yet relevant simplification we can cast the 
field-line Eqs. ([1]) in the Hamiltonian form 

dx dH{x,y,z) dy dH{x,y,z) , . 

dz dy ' dz dx 

by simplifying the magnetic field model to B(r) = Boz+6'Bj_{r±, z), where (5B_L(r_L, z) = 
VA||(r_L,z) X z is a static, transverse magnetic perturbation, A\\{r±,z) is the axial 
component of the vector potential, V = d/drj_, and r_|_ = {x, y) is the position vector in 
the perpendicular plane. Above, "time" is associated with the z coordinate, canonical 
momentum corresponds to the y coordinate, and x is the coordiante of a one and a half 
degrees of freedom system. With these settings the Hamiltonian is expressible in terms 
of the A\\ function as H{t±,z) = A\\{t±, z)/ Bq. Hereafter Bq = 1 for mathematical 
simplicity. The flow velocity is given by v_|_ = VH{r±, z) x z = 6^j_/Bq. If one 
introduces the Kubo number as i? oc v±/u!^±, one finds R oc {6B/Bq){^\\/^±), as 
expected, with t>_|_ ~ 6B/Bq ^ 1 a characteristic value of |v_l|, and u oc l/^y the 
"frequency" of turbulence variation due to the z dependence. 

3.2. Periodic orbits and separatrices 

To match with the numerical simulation model in Eq. ([2]), the vector potential is assumed 
to be given by a sum of plane waves, with the spectral components being odd functions 
of k^. At a given time z, consistently with the sign symmetry of the wave expansion, 
the separatrices of the A\\ field are given by 74||(r_L,z = const) = 0, while the closed 
equipotential contours are defined for A\\ 7^ as 

y = y{x, A|| = const, z = const). (5) 
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Figure 2. The Kolmogorov-Sinai entropy as a function of the Kubo number in the 
parameter range 1 < i? < 90, with the anisotropy level set to = 100- Dots with 

error bars: the numerically computed h values. Full line: a best fit by a logarithmic 
function (fit parameters in the legend). The inset shows the same data points, but 
the logarithmic fit (dash-dotted line) is done for the points up to i? = 40. Note the 
entropy slowing down for i? >~ 40. 



As usual, the area embraced by a closed contour is obtained as J = / ydx (this being the 
action of the Hamiltonian system). When the equipotentials in Eq. ([5]) are closed for all 
one speaks about a magnetic surface [3]. Magnetic surfaces are labeled by respective 
A\\ values, while / = /(^y) defines the integrals of motion. A field line wrapping 
around a magnetic surface experiences a periodic motion in the (x, y) plane with the 
frequency fi(/) = dH{I)/dI = \VA\\{r±,z)\/\VI{r^,z)\ ~ v±/i, where i ~ |V/(r_L,-2)| 
is the perimeter of the equipotential contour. In writing the first term in the sequence, 
use has been made of the standard definition of the action-angle variables (J, w) and 
of the equation for dw/dt = ^{I) in Hamiltonian dynamics [l3j. Assuming circular 
equipotentials one immediately gets I{r^,z) oc r^^, |V/(r_|_, 2;)| oc r±, and / oc i"^. 
Here, r± represents the typical "radius" of a closed contour, in the same sense as the 
"diameter" of a dust grain. In the case of strongly shaped equipotential contours, 
however, the latter scaling relations may not hold true. Of particular interest are the 
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equipotentials characterized by self-similar geometry in a broad range of spatial scales 
best described as fractals [H]. The significance of such equipotentials in low-frequency 
turbulence is addressed shortly. For fractal equipotential contours the simple scaling 
/ oc £^ is generalized to / oc i'^^'^^ (the so-called "area-perimeter relation" [15]) where 
dh is the fractal dimension of the equipotential line and we have omitted redundant 
dimensional constants for simplicity. The fractal dimension dh is defined by the scaling 
i oc r'^ (instead of £ oc r_|_ in non-fractal geometry). The dh values are generally 
fractional larger than 1. Eliminating i by means of the area-perimeter relation one 
obtains oc v^/I'^'^/'^. 

3. 3. Resonances and the percolation property 

In the basic theory of Hamiltonian chaos it is well established that the existence of 
resonances has an important impact onto the dynamics pT| 1^ 113]. In the model studied 
here, resonances occur between the wrapping motion of magnetic field lines around 
the magnetic surfaces and the magnetic field variation due to the axial dependence 
of the fluctuation field. In the presence of local instability the matching axial and 
wrapping frequencies or the matching harmonics of these favor departure from the exact 
periodic, energy conserving motion thus giving rise to transport phenomena in phase 
space. In the limit of very low frequencies, the conditions for a resonance [TT] require 
the corresponding excursion periods to diverge. Mathematically, this can be satisfied 
for fractal equipotential contours ^45j, characterized by diverging lengths due to the 
structure on many scales. 

Next, because of the sign symmetry of the fluctuation field in Eq. ([2]), the zero-set 
A\\{tj^,z = const) = contains a percolating equipotential line, which stretches the 
entire system [161 HZ]- The percolating line is the channel through which the turbulent 
diffusion penetrates to the large scales. Analyses of Ref. [17| have shown that the 
equipotential contours in the model turbulence field ([2]) lying in close proximity to the 
percolating line A\\ are indeed characterized by the fractal geometry in the limit 
when the perpendicular correlation length — ^ oo, and their fractal dimension was 
found to be ci/i ~ 4/3, in agreement with the results of Ref. [48]. This fractal geometry 
is a generic property of fields with a broad power-law energy density distribution 
conventionally associated with "turbulence" [35]. Thus, when the fluctuations are very 
slow, we expect the resonance conditions to be naturally satisfied in vicinity of the 
percolating line, provided that the number of waves in Eq. ([2]) is sufficiently large. In 
order words, fractality emerging from the percolation property of y4||(rx, z = const) = 0, 
together with the wide spectral extension in Eq. ([3]), guarantees the existence of the 
resonant orbits for a; — thus supporting the turbulent diffusion in the limit i? ^ oo. 

3.4- Nonlinearity parameter 

A net result from the above reasoning is that the resonance conditions in the limit of very 
low frequencies (large Kubo numbers) are satisfied for fractal equipotential contours, 
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which, as is aheady mentioned above, are found in close proximity to the percolating 
line. The resonance conditions alone do not, however, define the microscopic behaviors. 
To this end, two important conditions need to be specified. One is the size of the 
nonlinearity parameter [111 HH] 



A = L 



dn{i) 



dl 



(6) 



In order for phase space diffusion to actually occur, this parameter must be large, i.e., 
A ^ 1. In this range, as was first observed by Chirikov, the resonances will overlap 
p!3] and the dynamics can, after memory for the initial conditions has been lost, be 
described statistically in terms of a Fokker-Planck equation |36] . Topologically, the loss 
of memory is due to the mixing of trajectories in phase space. The mixing time is 
obtained from Zc — l/2ci;lnA, while the number of overlapping resonances is ordered as 
VA > 1. 



3. 5. Chaotic and pseudochaotic properties 

Next the relationship between mixing and randomness must be addressed, which can be 
nontrivial in some cases. Indeed, because of the many overlapping resonances, the 
dynamics will, in fact, be random after a certain transient time, which one would 
estimate as, say, two-three characteristic times it takes to pass from one resonance 
to the other. When the phase space accessible for random motions is wide enough 
(ideally, comparable with the entire phase space of the system), then the transition to 
random dynamics occurs on time scales of the order of Zc. This is a typical situation 
in chaotic dynamics. It should be understood, however, that this typical situation does 
not always take place. Situations exist, in which the available phase space is just too 
narrow to permit good mixing properties. For instance, in low-frequency turbulence, 
random motions as dictated by pertinent resonance conditions can be squeezed to a 
subset of phase space with fractal geometry [2^. It is this a situation that occurs in 
critical percolation [16]. As a consequence, the dynamics only partially acquire the 
chaotic property in that they show random behavior, while the mixing turns out to 
be very poor (sub-exponential). This anomalously slow mixing may lead to "strange" 
phenomena in phase space, for instance to stickiness of phase space trajectories [26], just 
to mention one. The implication is that the Lyapunov exponents may vanish despite of 
the inherently random character of the motion. This type of behavior has come to be 
known as pseudochaotic [26l H9l 150]. In pseudochaos one finds, then, that the transition 
to random dynamics occurs on time scales much shorter than Zc- A simple model 
realization of this is found in random walks on percolation systems [23]. Anticipating 
our final conclusion below, here we mention that the observed properties of saturation 
of the KS entropy are explained in terms of the Hamiltonian pseudochaos. 
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3. 6. Order of limits 

The other relevant condition that one needs to identify the microscopic pattern 
of behavior brings together the two building blocks of the model: fractality and 
randomness. Focusing on the fractal properties we somehow expect the fractal range to 
be wide enough, ideally oo. On the other hand, the slowness of the fluctuations 

implies that R oc v±/ujC,± —>■ oo for uj —>■ 0. For the purpose of formal orderings, because 
the two limiting procedures may not be commutative, we require C,± ^ oo diverge faster 
than R ^ oo. That is, given a finite fluctuation frequency, we set S> v±/u!^±. The 
implication is that the fluctuations should not be too slow for the actual (finite) size 
of the system, otherwise the dynamics retain a deterministic character. Equivalently, 
^±/v± ^ l/u oc ,^11 , showing that the diffusion time to the distance must be large 
compared to the period of the field. Under these conditions, which are naturally satisfied 
for the small enough v± ~ 6B/Bq ^ 1, the dynamics will be inherently random with 
A ^ 1 even for very slow fluctuation frequencies. 

3.7. Scaling argument 

Let us now rewrite Eq. ([6]) in a more insightful form. Substituting oc v±/I'^'^^'^ 

and performing the trivial differentiation over I leads to A{£) oc v±/uji, where the area- 
perimeter relation has been applied. For the longest equipotentials present (r^ ~ ^^), 
the perimeter scales as £ oc making it possible to evaluate A oc v±/uj^'^\ 

Remembering that the Kubo number R oc v±/uj^± it is found that 

R/A^^i'-\ (7) 

This is a very important relation because it reveals that the Kubo number R can be 
very large if the geometry is fractal 1, dh > 1). When dh 1, fractal properties 

are lost, while the R and A parameters turn out to be of the same order. 

3. 8. Kolmogorov-Sinai entropy — functional form 

We are now in position to obtain the KS entropy for the system under consideration. 
In fact, the entropy, /i, can be evaluated as the inverse mixing time, z^. (called by 
some [H [15] the Kolmogorov length L^^to emphasize its geometric meaning), yielding 
~ l/^c — 2ci; In A (Ref. [H] for the full discussion). Eliminating A by means of Eq. ([71) 
it is found that 

/i~2u;lni?-2cu(4-l)ln^^-2cjlnC, (8) 

where the last term in the sequence arises from the suppressed dimensional constants 
in the above scaling relations. Given the anisotropy level, Eq. ([HI) defines the KS 

entropy in function of the fluctuation strength or the Kubo number, R. It is understood 
that, by changing i?, one also changes the contour or set of the contours on which the 
resonances occur. Fractal geometries of the resonant contours may also be different for 
the different R values. In this sense one should think of the dh parameter in Eq. ([8[) as 
being function of R. 
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3.9. Chaotic case 

The next step is to observe that the resonance conditions cannot simultaneously be 
satisfied for the fractal and non-fractal equipotential contours because of the immense 
difference in their lengths. The actual regime is determined by the degree of anisotropy 
of the turbulence, with stronger anisotropy favoring the implication of fractal properties. 
Consider first the non-fractal case; that is, assume the anisotropy be kind of moderate. 
Setting ci/i = 1 in Eq. (J?!), one finds i? ~ A. Thus, in the non-fractal regime, a large 
i? 3> 1 implies a correspondingly large A 3> 1, and, hence, the behavior is chaotic |46j . 
This conclusion finds support in the numerical simulation results of Ref. fl9\, where 
Poincare sections of field-line motion have been analyzed (Figs. 2 and 3 in their work). 
For dh = 1, the KS entropy in Eq. ([8]) reduces to ^ 2ujlnR — 2ujlnC. This logarithmic 
dependence on R is consistent with the numerical fit in the inset of Figure 2 for not too 
large the R values. 

3.10. Nonlinear saturation and pseudochaos 

With the increasing anisotropy (i.e., Kubo number) the resonant motions gradually 
shift to fractal equipotential contours characterized hj dh > 1. This causes the h value 
to deviate from a purely logarithmic growth, since the second term in Eq. ([8]) is now 
non-zero. Thus, the behavior is actually slower than a logarithm, as can be seen by 
close inspection of Figure 2. Clearly, the entropy growth slowing down occurs as a 
consequence of fractality infiuencing the resonance conditions. 

Furthermore, there is an upper bound on the entropy growth, and that is where 
the "saturation" [28] of the entropy comes into play. When the turbulence becomes 
extremely anisotropic, i.e., — ^ oo, the resonance conditions are satisfied in a very 
small neighborhood of the separatrix, A\\ 0. In this small vicinity, the resonance 
structure is actually very dense. If one defines the characteristic frequency, u, one finds 
the width of the resonance layer to be 6uj oc y/uj (as compared to the distance between 
the resonances Au oc u). See Ref. |Tl]. Note that 6uj ^ Au for — 0, and that the 
width of the resonance layer shrinks to zero (as a/cJ), while the density of the resonances 
in the layer, defined as inverse distance Au, diverges as 1/u. 

In a sense, uj being close to zero implies that all the resonances have overlapped 
in a very narrow stochastic layer enveloping the separatrix. Let the number of these 
resonances be A^max ^ 1- Since the number of overlapping resonances is given by V^, 
it is noted that the A value has reached its maximum; that is, Amax — (A/'max)^- Hence, 
for a given (finite) ^y, the KS entropy saturates at /imax — 2ti;lnAmax oc (4/^||) InA/'max- 
When oo, the resonance spectrum collapses, while the entropy /imax — > 0. 

Under the conditions specified above, the dynamics will inherently be random 
despite of the vanishing /imax- This means, precisely, that the asymptotic {R — > oo, 
^max 0) dynamics are pseudochaotic, in the sense of Refs. [261 1231 SHI [50] . To this 
end, our final conclusion is: The deviation from a purely logarithmic behavior and the 
ultimate nonlinear saturation of the KS entropy in the percolation limit are explained 
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by the natural tendency of the considered system to evolve into a dynamical state 
characterized by the pseudochaotic property. 

3.11. Connection to fractional derivative equations 

The dynamics, being pseudochaotic in the percolation limit, needs to be addressed. 
Because of the formally diverging mixing time Zc — l//imax oo the Fokker-Planck 
equation in the form of Ref. [16] can be written only asymptotically for z ^ Zc- 
On time scales shorter than Zc this conventional form of the Fokker-Planck equation 
is invalidated. There is a growing belief that, in many ways, pseudochaotic systems 
are described by space-time self-similarity and by kinetic equations written in terms of 
generalized (fractional) derivatives [26j. We note in passing that a fractional extension 
of the Fokker-Planck equation has, for Hamiltonian systems, been proposed in early 
works in Refs. [521 EI]. 

More so, because there is no unambiguous definition of the notion of fractional 
derivative, one important issue is to define in which sense the fractional operator is 
understood in the problem under consideration. In this spirit, Milovanov [23] has 
suggested that the percolation transport in Hamiltonian systems is described by a 
fractional time Fokker-Planck equation with the fractional derivative in the Caputo 
sense (as opposed to the alternative definition in the Riemann-Liouville sense) [53llS]. 
This suggestion conforms with the results of Refs. [551 ES] , where one finds a discussion 
of subtleties of definition of fractional derivatives and their relevance in plasma physics 
applications. 

3.12. Summary of the discussion 

Our main conclusion from the above analysis is: The entropy growth slowing down is due 
to the concentration of the random properties on a substrate of phase space with fractal 
geometry in response to the increasing anisotropy of the turbulence. The asymptotic 
behavior in the limit R oo and the paradigmatic concept of "percolation" transport 
[T6] conform with the notion of Hamiltonian pseudochaos (random non-chaotic dynamics 
with zero Lyapunov exponents) |26j . 




4. Conclusions 

In the present work, we have addressed the Kolmogorov-Sinai (KS) entropy, h, in 
turbulent diffusion of magnetic field lines in the parameter range of large Kubo numbers, 
-R S> 1. For not too large the R values (in practice, 1 < i? < 30), we found numerically 
that the h function is well approximated by a logarithmic dependence /i ~ p In i? -|- g 
(here, p and q are numerical fitting parameters). Most importantly, with the increasing 
R, the entropy logarithmic growth slows down showing signatures of saturation. These 
behaviors are in marked contrast with the prediction of Ref. [E] that h oc R^^^ In R for 




i?> 1. 
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To comprehend this apparent discrepancy, we have proposed a model of percolation 
transport of magnetic field lines, basing our investigations on general principles of 
Hamiltonian dynamics. We found that the logarithmic counterpart of the entropy is 
naturally explained by the Hamiltonian chaotic dynamics of magnetic field lines in the 
sense of Refs. fTR HHl [13] . The deviation from the purely logarithmic behavior and the 
slowing down of the entropy is explained by the concentration of the chaotic properties 
on a subset of phase space with fractal geometry. This fractal geometry is a generic 
property of fields with a broad power-law energy density distribution conventionally 
associated with the notion of "turbulence" [35]. Finally, the asymptotic saturation of 
the entropy is explained by the collapse of the resonance spectrum of the system in the 
limit u ^ (i.e., oo). 

At this point, a subtlety has occurred involving the order of limits (i.e., the fractal 
limit ^± — oo versus R oo). The convention adopted in this work is that ^± — >■ oo 
diverges faster than i? — oo. Under this last condition, the dynamics retain the random 
character despite of the vanishing KS entropy. This behavior bears signatures enabling 
to associate it with the notion of pseudochaos (random non-chaotic dynamics with zero 
Lyapunov exponents) [26l [IHl [50] . 

The final result of this investigation is that the paradigmatic "percolation" 
transport in the limit R oo involves two essential ingredients: randomness and 
fractality. Random motions are found to pursue in the limit a; ^ supported by 
the fractal geometry of the substrate on which the turbulent transport concentrates. 
It is the interplay of these two ingredients that explains the observed properties of 
behavior of the KS entropy. These analyses shed new light on the origin of percolation 
transport by associating it with a pseudochaotic (rather than simply chaotic) behavior. 
A formal treatment of the problem in terms of random walk models and fractional 
kinetic equations is given in Ref. [23] . 
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